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CHARACTERISTICS OF THE OPTICAL RADIATION FROM KAUFMAN THRUSTERS 
by Nelson L Milder and James S„ Sovey 
Lewis Research Center 

SUMMARY 

The optical radiation from plasma discharges of electron-bombardment mercury- 
ion thrusters was investigated. The purpose of the investigation was to determine the 
applicability of spectroscopic diagnositc measurements to such thrusters. To accom- 
plish this task, a spectrograph was used to determine the atomic and ionic species in 
the discharge plasma of a 30-centimeter-diameter hollow cathode thruster and a 1. 5- 
meter-diameter thermionic cathode thruster. A 0. 5-meter Ebert mounted optical 
monochromator was used to study spectral line intensity of radiation emanating from 
the discharge of the 30-centimeter-diameter thruster. 

The results of these measurements indicated that the discharge was composed pri- 
marily of mercury atoms and singly charged ions. Excitation spectra of doubly charged 
mercury ions was also detected in both the 30-centimeter- and the 1. 5-meter-diameter 
thrusters. 

Diagnostics performed on the 30-centimeter hollow cathode thruster resulted in 
(1) a theoretical interpretation of the observed spectra in terms of a superposition of 
primary electrons upon a Maxwellian distribution and (2) the identification of two mer- 
cury atom lines at 365. 5 and 365.0 nanometers (3655 and 3650 A) that were suitable for 
application of line amplitude ratioing techniques. The spectroscopic measurements indi- 
cated that electron temperatures in the discharge plasma ranged from about 1. 2 to 
7. 5 electron volts depending on discharge voltage and radial location. Also, the primary 
electron (defined as the ratio of primary electron density to Maxwell density) varied 
from near zero at three-quarters radius to about 5 percent on the thruster axis. 

Spectral line amplitude measurements in the neutralizer region indicated a 
neutralizer-beam coupling voltage of about 20 to 22 volts, in qualitative agreement with 
probe measurements of this potential for neutralizers operating in spot mode. These 
results suggest possible use of spectroscopic measurements to thruster control applica- 
tions by correlating detectable varying optical effects with changing thruster performance. 

The percentage of doubly ionized mercury with respect to the atom density in the 
discharge could be measured directly and was found to be a maximum of about 2 percent 
on the thruster axis for a discharge voltage of 58 volts, dropping to about 1. 3 percent at 
half radius. At lower discharge voltages, the ratio of double ions to neutral atoms did 
not exceed about 0. 8 percent. 



INTRODUCTION 


The purpose of this investigation was to determine if spectroscopic diagnostic 
measurements could be made in mercury bombardment thrusters. The discharge 
plasma in the ion chamber of a mercury- electron bombardment thruster is formed 
primarily by the collision of monoenergetic and thermal electrons with mercury atoms. 
Phenomenological descriptions of the plasma have been given elsewhere (refs. 1 and 2). 
In general, the plasma electron densities in such thrusters are of the order of 10*^ to 
10 11 per cubic centimeter (refs. 3 and 4). Maxwellian electron temperatures generally 
do not exceed about 10 electron volts (ref. 4), and the neutral atom density is estimated 
to be of the order of 10 per cubic centimeter (ref. 3). The values of these quantities 
vary with thruster size and operating conditions. The ionization produced by electron 
bombardment is accompanied by collisional excitation of the atomic and ionic species 
in the discharge. Decay of these excitations produces radiation, as a visual inspection 
of an operating thruster immediately reveals. The present work describes investiga- 
tions of this radiation and extends the results of our preliminary measurements (ref. 5). 

There are several motivations for using spectroscopic techniques to study thruster 
plasmas. Previous diagnostics of the discharge have been based almost entirely on 
Langmuir probe measurements (refs. 1 and 2). Although these measurements have 
yielded plasma density, potential, and electron temperature profiles, there is always 
the concern that electric probes can introduce local disturbances that affect the param- 
eters being measured. Spectroscopy is a diagnostic method that does not disturb the 
plasma. 

In addition, there are properties of the thruster plasma that are particularly con- 
ducive to spectroscopic analysis. The low particle density in the discharge yields a 
plasma that is optically thin except for resonant states. This means that the radiation 
emanating from processes occurring deep within the plasma volume can be observed. 
Thrusters can be maintained at a steady operating condition for long times compared 
with those needed to obtain spectra. The operating point can be readily changed to 
another stable condition, so that the effect of varying thruster parameters on the radia- 
tive processes can be studied in detail. 

There are limitations, however. Data interpretation is generally not direct, re- 
quiring a knowledge of or assumptions on the nature of the physical processes in the 
discharge producing the radiation. Accurate measurements require detailed information 
on the optical properties of the radiating atom. Often such information is not available 
for heavy atoms such as mercury. Also, because the emitted light emanates from 
atomic sources located along the line- of- sight of the detecting system, local properties 
of the discharge plasma cannot easily be measured. Thus, one generally obtains infor- 
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mation on the average properties of the discharge. Even with these restrictions it was 
found that data could be interpreted in terms of a physical model to be described herein. 


EXPERIMENTAL APPARATUS AND PROCEDURE 
Thrusters 

The 30- centimeter- and the 1. 5- meter- diameter thrusters used in this study are 
shown schematically in figures 1 and 2. Detailed descriptions of the 30- centimeter 
thruster have been given in references 6 and 7. In this thruster a hollow cathode serves 
as the source of ionizing electrons. These electrons are contained in the discharge 
chamber by means of field- shaping permanent magnets. Propellant utilization ranged 
from about 80 to 95 percent. (The 30-cm thruster was operated by V. K. Rawlin and 
R. T. Bechtel of Lewis for this investigation. ) 

The 1. 5- meter (anode) diameter thruster (ref. 8) uses 10 oxide- coated thermionic 
emitters to produce ionizing electrons. A very weak magnetic field (of the order of 


Permanent magnet 

/ 



Figure 1. - Section view of a 30-centimeter-diameter thruster. 
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Figure 2. - Section view of 1. 5-meter-diameter thruster. 


r Vacuum tank 



Figure 3. - Experimental arrangement. 


_ A 

10 T) provided the means for suppressing electron radial diffusion to the anode. Pro- 
pellant utilization was nominally 65 percent. (This thruster was operated by 
S. Nakanishi of Lewis. ) 


Facility 

The thrusters were rim in the Lewis Research Center 7. 6- meter- diameter by 
18. 3-meter long vacuum facility described in reference 9. Their position in the vacuum 
tank is shown in figure 3. Facility constraints required that the optical detectors be 
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located at the opposite end of the facility from the thrusters. This meant that the dis- 
charge chamber was viewed through the ion exhaust beam. Radiation emanating from 
regions downstream of the discharge chamber would be due to neutralizing electron- 
atom excitations, electron-ion recombination, and possibly charge exchange. The down- 
stream radiation was measured and found to represent only a negligible fraction of the 
total measured radiation. 


Optical Measurements 

Two optical detection systems were used. A spectrograph using photographic plates 
was useful for general surveys of the thruster spectrum and for excited species identi- 
fication. This system is particularly useful for detecting very weak transitions. A 
monochromator with a photomultiplier detector was used for quantitative measurements 
of individual spectral lines. 

Spectrograph . - A plane grating spectrograph of a Czerny- Turner design with 
3/4-meter focal length was used. The spectrograph was equipped -with a cassette for 
holding 12.7- by 17. 8- centimeter photographic glass plates. The cassette could be 
positioned by means of a vertical rack, so that several spectra could be obtained on a 
single plate. A mechanical shutter at the entrance slit was used to obtain exposure 
times ranging from 1/200 second to several seconds. 

The range of wavelengths studied was from 230 to 670 nanometers (2300 to 6700 A), 
with a dispersion varying linearly from 2. 1 nanometers per millimeter (21 A/mm) at 
370 nanometers (3700 A) to 1.93 nanometers per millimeter (19.3 A/mm) at 630 nanom- 
eters (6300 A). The best resolution obtained at this range of dispersion was of the order 
of 0. 15 nanometers (1. 5 A). The high speed (f/6.3) of the instrument was especially 
useful in obtaining spectra of the 30- centimeter- diameter thruster. The long distance 
between this plasma source and spectrograph entrance slit resulted in many lines having 
rather low intensity. A typical densitometer trace of photographed spectral lines is 
shown in figure 4. 

Monochromator . - An Ebert mounted, plane grating monochromator with a 0. 5- 
meter focal length was used to study discrete spectral line amplitudes. The instrument 
dispersion was 3.2 nanometers per millimeter (32 A/mm). The plane grating was 

o 

rotated to cover the wavelength range from 230 to 800 nanometers (2300 to 8000 A). The 
recorded photomultiplier signal was corrected for the spectral response of the optical 
system. This response was determined using a radiance standard calibrated by the 
National Bureau of Standards. 

The monochromator was used to obtain radial profiles of the 30- centimeter- 
diameter thruster. The procedure consisted of setting the monochromator on the peak 
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Figure 4. - Densitometer trace of spectral lines. 
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Figure 5. - 30-Centimeter-diameter thruster image 
viewed from downstream of the extraction grids. 


of a given spectral line. A motor driven 25- micrometer horizontal slit was then used 
to scan the image on the vertical entrance slit of the monochromator as shown in fig- 
ure 5. In this manner the area of thruster cross section observed represented about 
0. 005 percent of the total thruster cross sectional area. For measurements of the 
thruster discharge chamber plasma, profiles were obtained along the radial direction 
away from the neutralizer location. 


SPECIES IDENTIFICATION AND OBSERVED LEVEL EXCITATIONS 

A tabulation of most of the observed Hg I and Hg II transitions are given in tables I 
and n, respectively. Also given are the energies above ground state of the upper levels 
of the transitions. Identification was based on compilations given in the literature 
(refs. 10 to 14). Energy level diagrams are given in figures 6 and 7. 

The mercury atom contains 80 electrons in the configuration 



TABLE X. - OBSERVED Hg I SPECTRA 


Series 

Transition 

Wavelength, 

Upper state 

Series 

Transition 

Wavelength, 

Upper state 


(ref. 15) 

nm 

level, 


(ref. 15) 

nm 

level, 



(ref. 13) 

eV 



(ref. 13) 

eV 

Singlet 

X S 

- 6 1 P 1 

491.60 

9.23 

Triplet 

3 S 

7 3 S x - 6 3 P q 

404.66 

7.73 



9 x s 0 - 6^ 

410.80 

9.73 



7 3 S x - 6 3 P 1 

435.84 

7.73 



lO^p - 6 1 P 1 

380.17 

9.46 



7 3 S x - 6 3 P 2 

546.07 

7.73 








8 S S 1 - 6 3 P„ 

334.15 

9. 18 


1 P 

a 1 ?, - 7 x s f , 

671.62 

9. 55 








1 u 





8 J s 1 - e^p,. 

275.28 

9. 18 



9 1 P 1 - 7 1 S„ 

623.44 

9.92 








1 0 





8% - 6^ 

289.36 

9.18 



KTP. - 7 x S n 

580.36 

10. 05 








1 u 





9 d s.. - 6 3 p 0 

292.54 

9. 70 



11 1 P 1 - 7 1 S Q 

554.98 

10. 15 












3 p 

9 3 P n - 7 3 S 1 

587.20 

9. 86 


X D 

6 X D 0 - 6 X P, 

579.06 

8.86 








2 1 





g^p, - 7 ,5 s 1 

585.94 

9. 86 



7 1 d„ - e 1 ?. 

434.75 

9. 57 








2 1 





9 6 P 0 - TS. 

582. 15 

9. 87 



8 1 d„ - a 1 ?. 

390.64 

9.90 








2 1 





lO^P, - 7 J S 1 

535. 40 

10. 02 



9 1 d„ - e 1 ?. 

370.42 

10. 05 








2 1 





11 3 P 2 - 7 3 S x 

512.05 

10.15 

Singlet- 

l,3 g 

7 1 S q - 6 3 P X 

407.78 

5. 15 



12 3 P 2 - 7 3 S x 

498. 06 

10.20 

triplet 


9 3 S 1 - 6 1 P 1 

414. 04 

9. 70 












3 D 

6 3 d 1 - 6 i P 1 

313. 15 

8. 85 


1 P 

8 X P X - 7 3 S x 

607.26 

4. 55 



6 3 D 2 - 6 3 P 1 

312.57 

8. 86 



9 X P X - 7 3 S x 

567.59 

9.92 



6 3 D x - 6 3 P 2 

366.29 

8.86 



10 1 P 1 - 7 3 S x 

531.67 

10. 05 



6 3 D 9 - 6 3 P, 

365.48 

8. 86 



U 1 P 1 - 7 S S X 

510.24 

10. 15 



6 3 D 3 - 6 3 P 2 

365. 01 

8. 86 








6 3 D. - 6 3 P„ 

296.73 

8. 86 


3 p 

6 3 P, - 6 1 S„ 

253.65 

4. 89 








1 0 





7 3 D 2 - 6 3 P 2 

302.35 

9. 56 


3 D 

6 3 D x - 6 1 P 1 

578.97 

8. 86 



7 3 D 2 - 6 3 P X 

265.20 

9. 56 



6 3 D 2 - 6 1 P 1 

576.96 

8. 86 



7 3 D 3 - 6 3 P 2 

302. 15 

9. 56 



7 3 Dj - 6 1 P 1 

433.92 

9. 56 



8 3 D x - 6 3 P 2 

280. 45 

9.90 



8 3 D, - 6 1 ?, 

390. 19 

9.90 



8 3 D 3 - 6 3 P 2 

280.35 

9.90 








9 3 D 3 - 6 3 P 2 

269. 88 

10.0 
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TABLE H. - OBSERVED Hg II SPECTRA 


Wavelength, 

Transition 

Upper state 

Wavelength, 

Transition 

Upper state 

nm 

(ref. 15) 

level, 

run 

(ref. 15) 

level, 

(ref. 14) 



eV 

(ref. 14) 



eV 

284.77 


fi 2 p 

11.9 

448. 75 




7 5 1/2 

6 P 3/2 




319. 10 
320.82 

92 °3/2 - 

7^P 
1 *1/2 

17.3 

460. 66 
466. 03 

ll2p 3/2- 
72 °3/2 - 

s25 1/2 

62p i/2 

17. 7 
15.5 




325.22 

326.41 

1q2d 5/2 - 

7 2 p 

1 j/2 

17.7 

470. 46 
476.22 

8 G 7/2, 9/2 ‘ 

8 2 g 7/ 2 - 

52f 7/2 

52f 5/2 

17.9 

17.9 




331.23 

D 5/2 ' 

2 P 

16.7 

479. 70 




P 3/2 

9 2s 1/2 - 

72p 3/2 

16. 4 

338.52 


485. 57 



345. 17 

512.84 

72 °5/2 - 

® 2p 3/2 

15. 4 




349.83 

82 °5/2 - 

2 P 

16.7 

520. 48 




P 3/2 




352.42 


62 °3/2 

16.5 

521.08 




6 F 5/2 ' 

72g 7/2,9/2 ' 

52f 7/2 

17.6 

353.26 



522.28 



354.94 

529, 40 

72g 7/2 - 

9 2f 5/2 

17. 6 




359.35 

62f 7/2- 

62 °5/2 

16.5 

542.52 

52f 5/2 - 

62 °3/2 

15.3 

360.58 

9 2d 3/2 - 

72p 3/2 

17.3 

559.54 

52F 5/2 - 

62 °5/2 

15.3 

377.63 

383.93 

9 2s 1/2- 

2p 3/2 

16.4 

585.19 

612.33 

62g 7/2 - 

2 -p 

F 5/2 

17.2 







398. 40 

R 2p 

62 °5/2 

7.51 

614.69 




6 P 3/2 

72p 3/2 - 

72s 1/2 

13.9 




614.95 





421.25 

1q2g 7/2, 9/2 ' 

52f 7/2 

18.2 

629. 13 
639. 49 

® 2g 7/2, 9/2 - 

6 2 g 7 /2 - 

5 2f 7/2 

52f 5/2 

17.2 
17. 2 

439.86 

82 °5/2 - 

72p 3/2 

16.7 

650. 15 







442.52 

82 °3/2 - 

72p 3/2 

16.7 

652. 11 

82s 1/2 - 

62p 3/2 

15. 1 

445.55 

9 G 7/2, 9/2 " 


18. 0 

664. 67 




5 F 7/2 
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Figure 6. - Energy level diagram for Hg I. 


Figure 7. - Energy level diagram for Hg II. 



The two electrons in the outer shell give the atom a helium- like structure. That is, the 
atomic levels consist of singlet and triplet series. Because of the larger atomic number 
for mercury, however, the radiation spectrum is more complicated than that of helium. 
This complication arises from the appreciable splitting of the triplet levels and from the 
allowed singlet- triplet transitions in mercury. 

The energy level diagram of figure 6 includes the significant levels observed in this 
study. This diagram was constructed from the data of table I with the aid of refer- 
ence 12. With the exception of the 185- nanometer resonant transition (broken line), the 
transitions given in the figure are representative of those observed and recorded in 
table I. The 185-nanometer line was not observed because transitions at wavelengths 
less than 230 nanometers were not within the range of the optical detection system. 

Because the Hg II ion has only one electron in the outer shell it has a doublet struc- 
ture for the excited states. All the ion levels observed in this study are given in figure 7. 
Levels for which the electronic configurations were unknown are represented by hori- 
zontal dashed lines. The transitions (shown by long-dash lines) were not observed but 
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are presented for completeness. The lowest excited state of Hg II is a 6 meta- 

stable level with an energy of 4. 5 electron volts above the ground state. 

In addition to singly charged mercury ion excitation, excitation of doubly charged 
ions (Hg HI) was also observed. Excitation of Hg HI was detected most strongly at 
479. 7 nanometers. The upper level for this transition is located at 15. 6 electron volts 
above the ground state for the doubly charged ion. With long exposure times, other 
weak Hg HI excitations were observed at 331. 2 nanometers at the upper state energy of 
16. 5 electron volts, 421. 67 nanometers, at 15.0 electron volts, 355.7 nanometers at 
16. 5 electron volts, and 521. 0 nanometers at 15. 2 electron volts. Examination of the 
1. 5-meter-diameter thruster spectrum also revealed a very weak excitation at 
396. 8 nanometers, which Schaffernicht (ref. 15) has identified as Hg IV. 

Impurities observed in the discharge are tabulated in table HI. In general, only the 
most persistent lines of the particular element were observed. Only iron lines were 
observed in the 30- centimeter- diameter thruster, but barium, chromium, and iron were 
found as faint lines in the 1. 5- meter- diameter thruster. The barium originated from 
the oxide cathodes used. Chromium and iron probably originated from the stainless 
steel thruster housing and anode in the 1. 5-meter thruster. The iron in the 30- 
centimeter thruster could have originated from the thruster housing and anode or from 
the mild-steel hollow cathode pole piece. If the impurity source were mild steel, this 
might explain the absence of chromium in the spectrum of this thruster. 

The present study indicated that the excitation of impurities in the discharge plasma 
was small. The total number and intensity contribution of the impurity lines to the 
radiation spectrum at the normal thruster operating discharge energy was substantially 
less than that of either the Hg I or Hg II contributions. This result suggests that the 
effect of these impurities on thruster performance was negligible. 
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TABLE HI. - IMPURITIES IN 1.5- METER- 


DIAMETER THRUSTER 


Wavelength, 

Impurity 


Wavelength, 

Impurity 

run 



nm 


350. 11 

Bal 


425. 43 

CrI 

358. 12 

Fel 


427. 17 

Fel 

a 371. 99 

Fel 


427. 48 

CrI 

a 373. 71 

Fel 


428.97 

CrI 

374.94 

Fel 


455. 40 

Bal 

385.94 

Fel 


553.56 

Bal 




614. 17 

Bal 


a Also observed in 30-cm-diam thruster. 


THRUSTER DIAGNOSTICS USING SPECTRAL LINE AMPLITUDES 

Theory 

The observed spectral amplitude is a measure of the emitted photons of energy 
corresponding to the wavelength of the spectral line. The problem is to relate the ob- 
served number of transitions at a given wavelength to the physical processes populating 
the energy level producing the photon emission (ref. 16). The following statements can 
be made on the physical nature of the discharge plasma within electron bombardment 
thrusters. 

_ q 

The discharge operates at low pressure (less than 10 torr) so that collisional 
broadening of spectral lines, collisional energy transfer between atoms, and three- body 
volume recombination would be expected to have a negligible effect on the observed line 
amplitudes (ref. 17). With exception of the resonant states, absorption should play a 
negligible role in excited state populations. 

The ion chamber plasma, although essentially a steady- state discharge, was not in 
thermal equilibrium. This was shown by presuming that such equilibrium did exist in 
the discharge plasma and then applying equilibrium theory and the Boltzmann-Saha equa- 
tions (refs. 18 and 19) to calculate electron temperatures from the observed line ampli- 
tude. Electron temperatures T g were calculated that were one to two orders of mag- 
nitude below expected values based on Langmuir probe measurements. These low T g 
values result from application of thermal equilibrium equations to nonequilibrium sus- 
tained plasmas (ref. 20). 

Langmuir probe measurements of such discharges (refs. 2, 4, and 21) indicate that 
the electron distribution consists of nearly monoenergetic, or primary, electrons at an 
energy determined by the discharge chamber potential difference and a low- temperature 
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Maxwellian distribution. In order to interpret the observed spectral line amplitudes in 
terms of the electron energy producing the emitted radiation, it is necessary to incor- 
porate a theoretical model of the discharge plasma. Because, in principle, both kinds 
of electrons can contribute to excitation and ionization of mercury atoms, a theoretical 
model of the discharge plasma should incorporate excitation and ionization by both 
groups . 

With these preliminary considerations, the rate of populating the j** 1 excited state 
due to collision excitation from the ground state is given by 

N j = N o [ N m S j< T e> + ViM cm ~ 3 sec_1 (1) 

J-T_ 

where S, (T ) is the Maxwell averaged excitation coefficient for the j n state and S.(E ) 
J e IP 

is the excitation coefficient averaged over the primary electron distribution. (All 

symbols are defined in appendix A. ) 

The number of j — k transitions is then given by the steady- state relation 

N jk = N j A jk cm ~ 3 sec_1 (2) 

where A^ is the relative transition probability for the j — k transition, 

A^ = Aj^/Aj Equation (2) can be rewritten with the aid of equation (1) such that 

N jk ' N o N m [l + ? VT e , E pH W < 3 > 


where 


<P 


N m 


X (T E ) - p 

X jk U e’ V “ o (r j 

jk T e 


( 5 ) 


and 


S jk = <Qjk V e> 

where the brackets denote averaging over the appropriate electron distribution. The 
quantities Q,, are related to the excitation cross sections Q. by the relation 
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Q jk - Q j A jk 


( 6 ) 


They are the optical excitation functions for producing the observed spectral line cor- 
responding the the j — k transition. The most recent measurements of optical excita- 
tion cross sections for mercury are those of Anderson, Lee, and Lin (ref. 22). A com- 
plete compilation of available mercury atom optical excitation functions can be found in 
reference 23. 

In order to obtain the electron temperature and primary electron fraction, it was 

necessary to form a ratio of two line amplitudes. Examination of the observed mercury 

spectrum resulted in the selection of spectral lines at 365 nanometers (3650 A) 

(6 3 D 3 - 6 3 P 2 ) and 365.5 nanometers (3655 A) (6 3 D 2 - 6 3 P 2 ). Population of the 6 3 D 3 

level from higher lying states by cascading is negligible in the electron- energy range of 

3 

interest (ref. 22), and there is no singlet mixing of this state. The 6 D 9 , however, 

1 " 
mixes strongly with the 6 D 2 to produce an optical excitation characteristic of the 

singlet state (ref. 22). These optical excitation functions are shown in figure 8. The 

calculated 365. 0- nanometer excitation cross section using Gryzinski’s theory (refs. 24 

and 25) is also given for purposes of comparison. The excitation cross sections for 

these two lines are thus sufficiently different in their dependence on electron energy 

that the line ratios vary significantly. Thus, the method described in references 20 

and 26 can be used to obtain electron temperature and primary electron fraction. 

From equation (3) the ratio of the number of 365. 5-nanometer transitions to 365.0- 

nanometer transitions (the ratio of spectral line amplitudes) is given by 


N 365. 5 
N 365. 0 


1 + yX 365 5 (T e ,E p ) 
} + ^ x 365. 0^ T e’ E pl 


S 365.5( T e) 
S 365. O^V 


250 xl 0' 19 



(a) For 6^03-6^2 transition. (b) For 6^D2-6^P 2 transition. 

Figure 8. - Optical excitation functions for Hgl (refs. 24 and 25 ). 
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The variations of the appropriate excitation coefficients with electron temperature or 
primary electron energy are given in figures 9 and 10. 

In order to use equation (7) to obtain primary electron fractions <p and electron 
temperatures T , it was necessary to choose a value for the primary electron energy 
E . Previous studies of the discharge of such thrusters (ref. 27) indicated that the 
primary electron energy can be related to the discharge chamber potential difference by 
the expression 


Ep = AVj - V p (8) 

where Vp is the plasma potential in the cathode- magnetic pole piece region. In the 
present study this potential Vp was taken to be 16 volts, based on electric probe meas- 
urements reported in reference 27. 

A qualitative estimate of the sensitivity of the transition number ratio 
■^365 5^365 0 *° P r i mar y electron energy is shown in figure 11. Here, the left side 
of equation (7) was plotted against E at two electron temperatures and three values of 
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Figure 11. - Theoretical transition number ratio against primary 
electron energy. 
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primary fraction. It was evident from such plots that the transition number ratio be- 
came less sensitive to changes in primary energy as the electron temperature was in- 
creased. Thus, an error in the selection of E would have the largest effect at the 
lower electron temperatures and high primary fraction. 


Atom Excitation 

The preceding theory was used to interpret the results of spectroscopic measure- 
ments obtained on a 30- centimeter- diameter hollow cathode bombardment thruster. 
Spectral line amplitudes were obtained for discharge chamber potential differences AVj 
of 33, 38, 43, and 58 volts. The discharge current was maintained at 8 amperes, and 
the thruster beam current was constant at 1. 5 amperes. At each discharge voltage the 
ratio of the number of 365. 5- to 365. 0- nanometer transitions was calculated at four 
radial positions. These ratios are shown graphically in figure 12. Here the measured 
ratio Nggg 5 /N 305 q is plotted at the four radial positions. Of particular interest is 
the minimum occurring at 22 electron volts at all radial positions except that furthest 
from the thruster axis. Such a minimum can be explained by the existence of nonran- 
domized primary electrons in the plasma discharge with the aid of figure 13 (see appen- 
dix B). 

The results of the analysis are summarized in figure 14. Here, the variation in the 
range of electron temperature with primary electron energy is presented for four radial 
positions. The boundaries of the crosshatched regions were determined by the range of 



Figure 12. - Measured transition number ratio against primary electron energy, 
where Ep ■ AVj - 16 eV. 


50 
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(c) Primary electron energy, 28 electron (d) Primary electron energy, 42 electron 

volts. volts. 


Figure 13. - Theoretical transition number ratio against electron temperature. 


primary electron fraction. The values of electron temperature obtained were in general 
agreement with Langmuir probe measurements of hollow cathode thruster discharges 
(ref. 4). The result that the largest values of primary electron fraction and electron 
temperature were obtained on the thruster axis, whereas the smallest values of these 
parameters were obtained at three-quarters radius, is consistent with recent measure- 
ments of a primary electron region centered about the thruster axis (ref. 28). 


Effects of Metastables 

For mercury the contribution to excited states by collisional excitation from meta- 
stable states can be shown to be less than the direct excitation from the ground state and 

o 

could thus be neglected. The ratio between collisional excitation of the 6 state from 
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(c) Normalized radial position, 0.5. 


(dl Normalized radial position, 0.75. 


Figure 14. - Electron temperature against primary electron energy. 


o 3 

the 6 Pq metastable level to the rate of collisionally populating the 6 D 3 from the 
ground state was calculated to be about 0. 005 over an electron temperature range of 4 to 
7 electron volts. The corresponding ratio for excitation from the 6^2 metastable 
state was about 0. 027 over the same electron temperature range. These ratios were 
even lower at electron temperatures below 4 electron volts. The calculations were 
based on the Gryzinski theory (refs. 24 and 25) and the assumption of a pure Maxwellian 
electron distribution. Experimentally, it was found that the primary electron density 
was generally less than 5 percent of the Maxwellian electron density. This fact, in con- 
junction with the fact that excitation coefficients for primary and Maxwellian electron 
excitation are of the same order of magnitude (figs. 9 and 10), implies that these esti- 
mates of the metastable level contribution to equation (1) are reasonable. Thus, exclu- 

3 

sion of the metastable contribution to the collisional excitation of the 6 C states was 
justified. 


Ion Excitation 

Table IV lists the electron energy required to produce various degrees of ionization 
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TABLE IV. - IONIZATION ENERGIES AND CROSS SECTIONS 


FOR MERCURY (REFS. 12, 23, AND 28) 


Ionization reaction 

Minimum energy required, 

Maximum cross section, 


eV 

cm ^ 

Hg I — Hg n 

10.43 

5. 4X10' 16 

Hg I - Hg m 

29.2 

5. 4X10' 17 

Hg I - Hg IV 

63.4 

3. OXIO' 18 

Hg n - Hg m 

18.75 

4. 2X10' 16 

Hg HI - Hg IV 

34.2 

6. Otx 10' 16 


of mercury by electron bombardment (ref. 12). Kupriyanov, and Latypov (ref. 29) have 
measured cross sections for the last two reactions of table IV, and several authors 
(ref. 23) have measured the cross section for production of multiply ionized mercury 
atoms. Optical excitation functions for mercury ions are essentially unavailable. In 
reference 29 the excitation function for the 284. 7- nanometer line of Hg II is reported. 
Here, however, excitation was the result of electron - neutral- atom collisions rather 
than electron- ion collisions. In the present work, calculated excitation cross sections 
based on the Gryzinski semiclassical method were used to obtain mercury ion excitation 
coefficients. 

For ion excitation the spectral amplitude at 398. 4 nanometers, corresponding to the 
2 2 2 

6 F 3 /2 “ 6s D 5/2 transition in n (see fig. 7), was used as the diagnostic line. It 

has been estimated (ref. 30) that contributions from cascading to the population of the 
o 

6 1 * 3/2 ^ eve ^ should be less than 10 percent. This line produced the strongest signal 

in the observed Hg H excitation spectrum. The calculated excitation cross section for 
2 

the 6 P 3/2 level is given in figure 15. Also shown is the calculated level excitation 
cross sections for the 5 djj /2 3/2 ®s^(J=4) Hg HI line at 479.7 nanometers. These cross 
sections are plotted as functions of electron energy. Excitation from the ground state 
of the ion (single ion for the 398. 4-nm line and double ion for the 479. 7-nm line) was 
assumed. From comparisons between such calculated cross sections for helium and 
experimentally measured helium optical excitation functions (ref. 25) the calculated 
cross sections for mercury are expected to be within an order of magnitude of actual 
cross sections. 

The variation of the ion to neutral fraction was determined as a function of primary 
electron energy at different radial locations in a 30- centimeter hollow cathode thruster. 
Using the method to obtain equation (7) measured amplitudes at 398. 4 nanometers for 
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Figure 15. - Theoretical level excitation cross sections for Hgll and Hg III. 


Hg II and 365. 0 nanometers for Hg I were ratioed to obtain the expression 


N 398. 4 1 + ^ X 365. 0^ E p* T e^ S 365.0 (T e^ 
N 365. 0 1 + ^ X 398. 4^ E p’ T e^ S 398. 4^ T e^ 


The quantity A 3gg 4 is the relative transition probability. Its value is unknown, but is 
expected to be substantially less than one. This is because the 6 2 P g ^ upper state of 
the 398. 4- nanometer transition is a resonant state and thus would most probably radiate 
to the ground state. The excitation coefficient coefficient ratios X x (E p , T g ) have been 
defined earlier and can be obtained from figure 16. Because A ggg ^ is constant, the 
right side of equation (9) yields the relative variation of N + /N q . The range of values of 
cp and T g used were obtained from figure 14. The relative variation of N + /N q is 
given in table V (third column) as a function of E p and radius. The calculations indi- 
cate that the ion fraction N + /N q did not vary appreciably with discharge voltage. 
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Figure 17. - Normalized ion and atom density variations with elec- 
tron energy. 30-Centimeter-diameter hollow-cathode thruster; 
ion beam current, 1.5 amperes. 


The relative variation of N + can be obtained from the relation 


N cc 

+ 



( 10 ) 


where it is assumed that N + ~ N m - Combining equations (10) and (9) enables us to 
obtain the relative variation of Nq. These variations are shown in figure 17, in which 
an average of the normalized ion and neutral atom densities are plotted as functions of 
primary electron energy. Maximum densities occurred at the lowest energy (17 eV), 
and all the data were normalized to unity at this energy. The marked decrease in both 
ion and atom densities can probably be attributed to the decrease in propellant flow re- 
quired to maintain the constant beam current (1. 5 A) with increased discharge chamber 
potential difference. This fact suggests the possibility of relating propellant flow 
changes to changes in atom and ion spectral amplitudes. Further studies at varying 
thruster operating conditions would be required, however, to determine whether spec- 
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troscopic techniques could be used to monitor propellant flow variations. 

Excitation of doubly charged mercury ions at 479. 7 nanometers were also measured 
in the 30- centimeter- diameter thruster. An expression analogous to equation (9) (with 
479. 7 nm substituted for 398. 4 nm) can be used to obtain N ++ /N q . The appropriate 
excitation coefficients are given in figure 18. In this case, A^g ^ = 1, because the 
upper state of this transition has only a single path for radiative decay (ref. 31). Thus 
the spectral amplitude ratio gives the fraction of N in the discharge directly. These 
fractions are given in the last column of table V. The accuracy of these measurements 
depends markedly on the accuracy of the excitation coefficients used. Because, to the 
authors' knowledge, no measured excitation functions for Hg m are available, it is 
impossible at present to assess the accuracy of these calculations. However, the 
N, /N ratios of table V were found to be consistent with measurements of double ion 

++ O 

currents obtained in the ion beam exhaust (ref. 32). 

It should be noted that excitation of Hg IH at an electron energy of 17 electron volts 

was observed at all radial positions, although quite weakly. The existence of Hg in at 

such low electron energies is to be expected because the threshold for ionizing Hg H to 
2 

Hg HI from the 6 metastable level of Hg II is about 15 electron volts. In addition, 

the "tail" of the low- temperature Maxwellian distribution of electrons can produce both 
ionization of Hg II and excitation of Hg HI. 


A POSSIBLE MEASUREMENT OF NEUTRALIZER-BEAM COUPLING VOLTAGE 

The orthogonal slit arrangement (fig. 5) was used to obtain the N 3 g 5 5^365 0 
ratio in the region of the neutralizer. The results are shown in figure 19. The center 
of the discharge chamber was located at r = 0. The N 3 gg 5^365 0 transition ratio 
in the neutralizer region was of the order of 0. 14 to 0. 18. These values were substan- 
tially greater than those obtained from the discharge region (fig. 12). Using the primary 
plus Maxwellian electron energy distribution theory described by figure 13, such high 
values suggest a large primary electron concentration and a very low electron tempera- 
ture (of the order of 1 eV). Hence, excitation by Maxwellian electrons can be neglected. 
A primary electron energy in the neutralizer region of about 20 to 22 electron volts was 
thus obtained. This result is in good agreement with potential probe measurements of 
the neutralizer- beam coupling voltage for neutralizers operating in spot mode (ref. 33). 
The Nggg 5/N305 q ratio dropped off rapidly toward the thruster side of the neutralizer, 
suggesting randomization of neutralizer electrons. The intensity ratios in the neutral- 
izer region appeared to be relatively insensitive to varying thruster parameters at a con- 
stant beam current of 1. 5 amperes. 
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TABLE VI. - TRANSITIONS FROM 
THE 7 3 S 1 LEVEL IN MERCURY 


.12 






Discharge 
voltage, 
A V If 

Emission 

current, 

J E . 

Wavelength, 

nm 

Transition 

0 

0 

0 0 

V 

O 32.5 

A 58.2 

A 

10 

7.9 

404.7 

435.8 
546. 1 

7 3 S 1 - 6 3 P q 
7 3 S x - 6 3 P 1 
7 3 S x - 6 3 P 2 


Thruster Neutralizer 


region 

1 


i region ^ 

1 

1 1 

10 


20 

30 40 


Radial position, cm 


Figure 19. - Measurement of neutralizer transition number ratio for 
30-centimeter-diameter thruster. Ion beam current, 1.5 amperes. 


RELATIVE AMPLITUDES OF THE 7 3 S X - 6 3 P 0) 12 TRANSITIONS 

Table VI lists the wavelengths of the so-called visible triplet of the mercury spec- 
trum. Because these transitions originate from the same state (7 S^), the relative 
spectral line amplitudes are independent of the level excitation cross section. The 
transition ratios 7^435 8 : ^546 1 are * hus a direct measure of the relative trans- 

ition probabilities for these lines. Several authors have reported measured values for 
these ratios, most of which have been tabulated in reference 34. In table VII, this tab- 
ulation is repeated, including the measurements of Anderson (ref. 22) and the present 
authors. In reference 22 the transition probabilities were calculated to obtain the theo- 
retical values for these ratios given in the table. 

The results of the present work were based on seven measurements made along 
different radial locations in the 3 0- centimeter- diameter hollow- cathode thrusters. The 
resulting relative transition probabilities are seen to be in general agreement with pre- 
vious work and in good agreement with the reported theoretical values. This fact pro- 
vides additional verification of the measurements obtained in this study. 
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TABLE VH. - RELATIVE TRANSITION RATIOS 


FOR THE VISIBLE TRIPLET IN Hg 


Source 

4047 

4358 

5461 

Hanle and Schaffernicht 

47 

100 

82 

(ref. 35) 




Theime (ref. 36) 

43 

100 

77 

Lebedeva and Fabrikant 

43 

100 

86 

(ref. 37) 




Schouten (ref. 38) 

37 

100 

134 

Jongerius (ref. 39) 

37 

100 

125 

Frisch and Klucharyov 

37 

100 

116 

(ref. 34) 
Anderson, et al. 
(ref. 22): 




Theory 

46 

100 

70 

Experiment 

36±3 

100 

96±13 

Present work 

38±4 

100 

71±14 


CONCLUDING REMARKS 

The excitation spectrum of mercury atoms and ions produced by electron bombard- 
ment thrusters has been determined. The discrete spectrum has been categorized for 
wavelengths between 230 and 700 nanometers (2300 and 7000 A). Excitation of singly, 
doubly, and possibly triply ionized mercury has been detected. 

Spectral line amplitude measurements indicated that the discharge was a nonequi- 
librium, sustained plasma. The energy distribution of electrons in the discharge of 
hollow -cathode thrusters can be represented by a distribution of monoenergetic elec- 
trons superimposed on a Maxwellian electron distribution. Ratios of the spectral line 
amplitudes at 365. 5 and 365.0 nanometers (3655 and 3650 A) were formed to obtain 
electron temperatures and ratios of primary to Maxwellian electron densities at four 
primary electron energies and different radial locations. It was found that the electron 
temperature ranged from about 1. 2 to 7. 5 electron volts and that the primary to Max- 
wellian electron density ratio ranged from zero to about 0.05. These values depended 
on primary electron energy and radial location in the thruster. 

Normalized ion and neutral densities decreased with increasing discharge voltage at 
constant emission current and beam current. It was presumed that this behavior was 
attributable to a decreasing propellant flow with increasing voltage. This decrease in 
flow was required in order to maintain a constant beam current. Spectroscopic meas- 
urements in the region of the neutralizer operating in spot mode indicated a coupling 
voltage of the order of 20 to 22 volts. 

The fraction of doubly charged mercury ions did not exceed about 2 percent of the 
atom density in the discharge plasma at a discharge voltage of 58 volts. At lower dis- 
charge voltages this fraction did not exceed about 0. 8 percent. 

One of the primary purposes of this investigation was to determine the extent to 
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which the optical radiation of electron bombardment thrusters could be used to study the 
discharges of such thrusters. The high sensitivity of the radiation output to changes in 
thruster parameters suggests potential engineering applications. For example, it may 
be possible to use the light output from the neutralizer region to control the neutralizer 
gas flow. It was concluded that this study demonstrated the applicability of spectros- 
copy as a diagnostic tool for thruster investigations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 19, 1971, 

113-26. 
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APPENDIX A 


SYMBOLS 


A ik 

transition probability for 

_ i 

n p 

primary electron density, cm 


j — k transition, sec 

n a 

number of transitions at wave- 

A j,tot 

total transition probability 


length X (spectral line am- 

for state j, sec - 


plitudes) 

A ik 

relative transition proba- 

N +,~ 

density of singly (and doubly) 

_ 3 

bility, A jk /A jjtot 


ionized atoms, cm 

a a 

relative transition proba- 

p 

atom orbital angular momentum, 


bility at wavelength X 


equal to 1 

D 

atom orbital angular momen- 

p 

electron orbital angular momen- 


tum, equal to 2 


tum, equal to 1 

d 

electron orbital angular 

Qj 

excitation cross section, cm^ 


momentum, equal to 2 


optical excitation function for 

E p 

primary electron energy 

2 

j — k transition, m 


(eq. (8)), eV 

r 

radial position relative to 

f 

electron orbital angular mo- 


thruster axis 


mentum, equal to 3 

S 

atom orbital angular momentum, 

Hgi,n,m 

mercury atom, single ion 


equal to 0 


and double ion, respec- 

s 

electron orbital angular momen- 


tively 


tum, equal to 0 

j 

total atom angular momen- 

S,(E ) 

primary electron excitation coef- 


tum 

J Jr 

ficient at energy E to state 

O 1 P 

J B 

ion beam current 


j, cur sec - 

N i 

rate of collisionally popu- 

1L 

lating j state from the 

Sj(T e ) 

Maxwell electron excitation 
coefficient at temperature T 


ground state 


3 1 ^ 

to state j, cur sec 

V 

number of j — k transi- 
-3 - 1 

S ik 

optical excitation coefficient for 

tions, cm sec 


j — k transition 

N m 

Maxwell electron density, 

q 

S A 

excitation coefficient at wave- 


cm 


length X 

N o 

-3 

neutral atom density, cm 

T e 

electron temperature, eV 
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Vp cathode- pole piece plasma poten- 

tial, v 

AVj discharge chamber potential dif- 
ference, v 

v g electron speed, cm sec” ^ 


W T e> 

X A<V T e> 


optical excitation coeffi- 
cient ratio, 


VVVV 

excitation coefficient ratio 
at wavelength X 


<P 


primary electron fraction, 

V N m 
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APPENDIX B 


METHOD FOR INTERPRETING THE OBSERVED ATOM EXCITATIONS 

In figure 13, equation (7) is plotted against electron temperature for the four pri- 
mary electron energies studied. Several curves corresponding to different values of 
primary electron fraction <p , ranging from 0 to 0. 1 are presented. It should be noted 
that curves of cp > 0 show a minimum. The intersection of the measured 
^365 5/N365 0 ra ^° wit* 1 these theoretical curves was used to obtain the ranges of 
electron temperature and primary electron fraction corresponding to a given value of 
primary electron energy E , making use of the following assumptions: 

r 

(1) Increasing the primary electron energy increased the available discharge power 
per atom and hence increased the electron temperature. Thus the electron temperature 
obtained from figure 13(d) would be greater than the temperature obtained from fig- 
ure 13(a). 

(2) Because an increase in primary electron energy increased the mean relaxation 
time for primary electrons, it was assumed that the primary electron fraction did not 
decrease with increasing primary electron energy. The average primary fraction ob- 
tained from figure 13(d) would thus be greater than the average fraction obtained from 
figure 13(a). 

(3) For values of primary electron energy E p less than 22 electron volts, the 
measured Nggg 5/^355 q ratio decreased with increasing electron energy. Also, the 
theoretical Nggg 5/ N 365 0 ratios figure 13 decreased with increasing electron tem- 
perature for temperatures below about 4 electron volts. Thus, for E^ < 22 electron 
volts, points of intersection of the measured Nggg 5/^355 g line to the right of the 
minimum of the theoretical curves of figure 13 were neglected. Similarly, for E^ > 22 
electron volts, only points of intersection to the right of the minima were retained. 

The allowable range of values for primary electron fraction at different radial loca- 
tions was obtained using condition (2) and the following considerations. At a given radial 
location (e.g. , r = 0) and primary electron energy equal to 22 electron volts, the meas- 
ured N^gg 5/^305 0 ratio (equal to 0. 065 at r = 0) was drawn tangent to the minimum 
of a particular cp- curve of figure 13(b). (For r - 0, the measured ratio at 0. 065 was 
tangent to the primary electron fraction, cp = 0. 03 curve at its minimum. ) The result- 
ing value of <p was the upper bound of the cp (and thus T ) range associated with 
Ep = 17 and 22 electron volts. The lower bound was obtained from the intersection of 
the measured Nggg 5/N305 q ratio at E^ = 17 electron volts with portions of the 97- 
curves of figure 13(a) to the left of the minima and from the fact that electron tempera- 
ture was presumed to exceed about 1 electron volt at all thruster conditions. This 
assumption was based on the fact that temperatures lower than 1 electron volt charac- 
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terize a thermal equilibrium plasma. From our earlier discussion, it was argued that 
the discharge chamber plasma was not in thermal equilibrium. Thus, for the radial 
position r = 0, where Nggg g/Nggg 0 = 0*074 E p = ^ electron volts, the lower 
bound on the primary electron fraction cp (from fig. 13(a)) was about 0. 01. 

For Ep > 22 electron volts, the procedure used to obtain the appropriate ranges 
of electron temperature and primary electron fraction was as follows: The maximum 
primary electron fraction was obtained from the intersection of the measured 
■*^365 5 /N 365 0 ra *i° with the curves of figure 13(d). For the radial position r = 0, 
the largest value of cp intersected by Nggg 5 /N 365 0 ec l ua i *° 0* 089 was about 0. 05. 
Thus, for primary energies less than 42 electron volts, intersection of measured 
N 365 5 /N 365 0 va -i ues with curves of cp > 0. 05 were disregarded, in accordance with 
condition (2). 

The minimum value of <p for E > 22 electron volts was the same as the minimum 
cp value for Ep < 22 electron volts, again, in accordance with condition (2). Thus, 
for r = 0 and Ep > 22. electron volts, the range of <p was found to be 
0. 01 s (p <. 0.05. For E £ 22 electron volts, the range was 0. 01 s cp < 0.03. 

Using the <p ranges and the appropriate curves of figure 13, the range of electron tem- 
perature for different values of E p could be determined. 
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